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The Quark Propagator in Momentum Space 
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The quark propagator is calculated in the Landau gauge at /3 — 6.0. A method for removing the dominant, tree- 
level lattice artefacts is presented, enabling a calculation of the momentum-dependent dynamical quark mass. 
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1. Introduction 

The quark propagator is one of the fundamental 
quantities in QCD. By studying the mass func- 
tion, which is the scalar part of the quark prop- 
agator, we can gain insight into the mechanisms 
of chiral symmetry breaking. The momentum de- 
pendence of the quark propagator is also used ex- 
tensively as input in Dyson-Schwinger equations 
for hadronic matrix elements. A lattice study of 
the quark propagator may enable us to check the 
validity of the models used in these calulations. 

2. C(a)-improved quark propagator 

All 0(a) errors in the fermion action can be re- 
moved by adding terms to the Lagrangian 



C(x) 



b g am 



+ ciaipp 2 ip + czamippip (1) 



improved propagator is given by 

S(x, y) = (i>'{xj$\y)) - L(x)So(x, y)R(y) (3) 

where So(x,y) = (ip(x)ip(y))- Since the propa- 
gator 5*o is defined as the inverse of the fermion 
matrix M(x) = p(x) + mo + O(a), we can use 
this to obtain another, simpler expression for the 
improved propagator: 

S(x, y) = (l + 2(b q + c q )am)S (x, y) 

- 2ac q S(x ~y) + 0{a 2 ) (4) 

With this in mind, we define the tree-level im- 
proved propagator Si(x, y) as 

Si(x-y) = (1 + ma)S (x-y) - -5{x-y) (5) 

We will also introduce the tree-level 'rotated' 
propagator Sr(x, y) as the special case of the im- 
proved propagator in (|^) , 

nun 

S R (x, y) = (l + —)L'(x)S (x, y)R'(y) (6) 
where L' = 1 — op /4 and R' = 1 + op /A. 



m should here be taken to be the subtracted bare 
mass m = m,Q — m c . The b g and b m terms corre- 
spond to a rescaling of the coupling constant and 
the mass respectively. The two last terms can be 
eliminated by a field transformation Q , 

ip -> ip' = (1 + b q am)(l - c q ap)ip = Lip 

= (l + b q am)Jj(l + c q a*p) = ^R (2) 



The tree level improved action after the transfer 
mation (0) has c sw 



— l,b q — j and c q = \. The 



3. Analysis 

The momentum space quark propagator in a par- 
ticular gauge is given by S(p) = J2 X g~ wx S(x, 0). 
We introduce the following 'lattice momenta' 
= isin(p M a) and = |sin(p^a/2) which 
differ by 



a 2 Ak 2 



0(a 4 



(7) 
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In the continuum, the quark propagator has the 
following general form, 



S(p) 



z c (p) 



i 0A c (p) +B c (p) i i> + M c (p) 



(8) 



We expect the lattice quark propagator to have a 
similar form, but with ji replacing p 1 : 



S(p) 



Zip) 



i$ + M (p) 



(9) 



The dimensionless Wilson fermion propagator at 
tree level is 



-i fta + mo<x + ifc 2 a 2 



(10) 



k 2 a 2 + (mod + ifc 2 a 2 ^j 
The tree level 'improved' propagator is given by 



S\ 0) (p) = (l + m a)S^(p)-- (11) 



Mi 



If we write 

we find 

Z(°>(p) 



i }ia + m a + aAM^ (p) 
Z(°)(p) 



fc 2 a 2 (l + m a) 2 + S 2 j 



(12) 



(13) 



(1 + moa)D 
aAM(Q)( P ) = TO ° a2 -f Afc2 + ^ 4 / 4 (14) 



1 + moa 



where 



L> = Fa 2 + ( 77i a + 2 fc2fl2 



m 2 a 2 a 4 Afc 2 a 4 /c 4 



S 



(15) 
(16) 



The corresponding tree level expressions for Sr 
are considerably more complicated than those for 
Si, and we will not reproduce them here. 

We are here primarily interested in studying the 
deviation of the quark propagator from its tree 
level value. Since QCD is asymptotically free, 
we should at sufficiently high momentum values 
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Figure 1. Z r (p) (left) and M r (p) (right) for all 
three propagators, from 20 configurations. 



get S(p) — > S(°'(p) up to logarithms. We will 
attempt to separate out the tree level behaviour 
by writing 



ia ft + aM r (pa) + aAM<°) (pa) 
Z r (pa)Z(°)(pa) 



(17) 



where mo is replaced by m in the expressions for 
and AM' '. Asymptotically, we expect that 
Z r (pa) — > 1 and M r (pa) — > m up to logarithms. 

4. Results 

The quark propagator is calculated at j3 = 6.0 on 
a 16 3 x 48 lattice, using the mean-field improved 
value c sw =XA79. The configurations were fixed to 
Landau gauge with an accuracy of 9 m&K — 10~ 12 . 
At k = 0.137, corresponding to ma = 0.0603, we 
have generated both So and Sr. At k = 0.1381, 
corresponding to ma = 0.031, only iSo was gener- 
ated. Si is easily constructed from So- 

When we calculate Z and M without factoring 
out the tree-level behaviour, it becomes clear 
that both improved propagators (Si and Sr) 
are completely dominated by the unphysical tree 
level behaviour at high momenta. In particular, 
B = ZM computed from Si becomes large and 
negative, approaching a value of aB = —3. Only 
in the infrared, below pa < 0.8, might we be able 
to extract physically significant information. 

Fig. [l] shows Z r and M r a as functions of pa for 
both Sq, Si and Sr. In order to ease the com- 
parison and remove residual anisotropy, we have 
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Figure 2. M r (p) from S7, for k = 0.137 (circles) 
and k = 0.381 (crosses). The increase in M r (p) 
for pa > 1 is an indication of the difficulty of accu- 
rately subtracting off the tree level mass function. 
The data shown are from 50 configurations. 



selected momenta lying within one unit of spatial 
momentum of the 4-dimensional diagonal. While 
all three show a dramatic improvement on the un- 
subtracted data, the large momentum behaviour 
of the mass function from Si is poor as a con- 
sequence of the pathological behaviour of Sf^ (p) 
at high momenta. This is due to the cancellation 
of large terms in the subtraction for the mass. It 
is therefore desirable to use the definition Sr for 
the improved propagator .[] 

Below pa ~ 0.9, the values for M r agree within 
errors for the two versions of the improved prop- 
agator. In particular, the value for the infrared 
mass M r (0) comes out the same. In contrast, So 
yields a mass which is 3-4er higher. 

Fig. H shows M r (p) calculated from Sj for both 
quark masses. We see that the infrared mass 
changes only slightly as the bare quark mass 
is halved, pointing to a dynamically generated 
quark mass of (300 ± 30)MeV in the chiral limit. 
Z r turns out not to depend on the quark mass. 



5. Conclusion and further work 

We have used two different definitions of the 0{a) 
improved quark propagator. We make use of 
asymptotic freedom to factor out the tree level 
behaviour, replacing it with the 'continuum' tree 
level behaviour Z(p) = l,M(p) — m. This tree 
level subtraction dramatically improves the data. 
The relatively poor behaviour of the tree level 
subtracted Sj can be put down to the large tree 
level finite-a effects which require fine tuning to 
subtract off correctly. 

For pa < 1, we see that M r (p) falls off with p as 
expected. The values obtained from Sr and from 
Si are consistent, while those for the unimproved 
propagator Sq differ significantly. We find that 
M r (0) approaches a value of 300 ± 30 MeV in the 
chiral limit. 

We also find a significant dip in the value for Z r (p) 
at low momenta. This is more pronounced for the 
improved propagators than for the unimproved 
one. The next step will be a quantitative study 
of the functional form for Z r and M r 

Finite volume effects have not been studied, but 
there is no sign of any anisotropy at low p, indicat- 
ing that finite volume effects are small. Repeating 
these calculations at a different lattice spacing is 
also essential especially to get reliable results for 
the quark mass. 
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1 A similar conclusion was reached by C. Pittori 



4. C. Pittori, iep-lat/9909086 



